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Abstract: On directed Barabasi-Albert networks with two and seven neigh- 
bours selected by each added site, the Ising model with spin S = 1/2 was seen not 
to show a spontaneous magnetisation. Instead, the decay time for flipping of the 
magnetisation followed an Arrhenius law for Metropolis and Glauber algorithms, 
but for Wolff cluster flipping the magnetisation decayed exponentially with time. 
On these networks the Ising model spin S = 1 is now studied through Monte 
Carlo simulations. However, in this model, the order-disorder phase transition is 
well defined in this system. We have obtained a first-order phase transition for 
values of connectivity m = 2 and m = 7 of the directed Barabasi-Albert network. 
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Introduction 

Sumour and Shabat [HI2] investigated Ising models with spin S = 1/2 on 
directed Barabasi-Albert networks jSj with the usual Glauber dynamics. No 
spontaneous magnetisation was found (and we now confirmed this effect), in 
contrast to the case of undirected Barabasi-Albert networks |3J |5J |B] where a 
spontaneous magnetisation was found lower a critical temperature which in- 
creases logarithmically with system size. More recently, Lima and Stauffer [7] 
simulated directed square, cubic and hypercubic lattices in two to five dimen- 
sions with heat bath dynamics in order to separate the network effects form 
the effects of directedness. They also compared different spin flip algorithms, 
including cluster flips [H], for Ising-Barabasi- Albert networks. They found a 
freezing-in of the magnetisation similar to PjJ |2] , following an Arrhenius law 
at least in low dimensions. This lack of a spontaneous magnetisation (in the 
usual sense) is consistent with the fact that if on a directed lattice a spin Sj 
influences spin Si, then spin Si in turn does not influence Sj, and there may 
be no well-defined total energy. Thus, they show that for the same scale-free 
networks, different algorithms give different results. Now we study the Ising 
model for spin S = 1 on directed Barabasi-Albert network and different from 
the Ising model for spin S = 1/2, the order-disorder phase transition of or- 
der parameter is well defined in this system. We have obtained a first-order 



1 



phase transition 
Barabasi- Albert 



for values of connectivity m = 2 and m = 7 of the directed 
network. 



spin-1 , m=2 




temperature 

Figure 1: Plot of spontaneous magnetization versus temperature for various 
network sizes. 



Model and Simulation 

We consider the Ising model with spin S — 1, on directed Barabasi- Albert 
Networks, defined by a set of spins variables S taking the values —1,0 and 
+1, situated on every site of a directed Barabasi- Albert Networks with N 
sites. 

The Ising interation energy is given by 

E = -jj2J2^ s k (i) 

i k 
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Figure 2: Plot of the Binder parameter Bi versus K for severals systems 
sizes (N = 250, 500, 1000, 2000, 4000, 8000 and 16000). In the Part (a) 
m = 2 and Part (b) m = 7. 

where fc-sum runs over all nearest neighbors of Si. In this network, each 
new site added to the network selects m already existing sites as neighbours 
influencing it; the newly added spin does not influence these neighbours. To 
study the critical behavior of the model we define the variable e = E/N 
and m = X^Li Si/N . ^From variable of the energy measurements we can 
compute, the average energy and specific heat and energetic fourth-order 
parameter, 

u{K) — [< E >} av /N, (2) 
C(K) = K 2 N[< e 2 > - < e >% v , (3) 

Bi(K) = [1 " a^^U (4) 

where K = J/ksT, with J = 1, and A;# is the Boltzmann constant. Similarly, 
we can derive from the magnetization measurements the average magnetiza- 
tion, the susceptibility, and the magnetic cumulants, 

m(K) = [< |m| >} av , (5) 

X (K) = KN[< m 2 >-< \m\ > 2 } av , (6) 

3 < m > z 
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Figure 3: Plot of the Binder parameter Bi versus 1/N for m = 2 (circle) 
and m = 7 (square), and severals systems sizes (N = 1000, 2000, 4000, 8000 
and 16000). 

where < ... > stands for a thermodynamics average and [...] a „ square brackets 
for a averages over the 20 realizations. 

In the order to verify the order of the transition this model, we apply 
finite-size scaling (FSS) jU]. Initially we search for the minima of energetic 
fourth-order cumulant in eq. (4). This quantity gives a qualitative as well as 
a quantitative description of the order the transition [HI]- It is known ^l] 
that this parameter takes a minimun value -Bj, m m at the effective transition 
temperature T C (N). One can show ^2] that for a second-order transition 
limAT^oo (2/3 — Bi jTn i n ) = 0, even at T c , while at a first-order transition the 
same limit measures the latent heat \e+ — e_|: 
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intercept=0.532(1) 
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Figure 4: Plot of the Binder parameter 2/3 — Bi(K) versus 1/N for m = 2 
(circle) and m = 7 (square), and severals systems sizes (N = 1000, 2000, 
4000, 8000 and 16000). 



^(2/3 - Bi ,_) = L^±^ns±±±r (8) 

A more quantitative analysis can be carried out throught the FSS of the 
specific heat C max , the susceptibility maxima Xmax and the minima of the 
Binder parameter Bi jTn i n . If the hypothesis of a first-order phase transition is 
correct, we should then expect, for large systems sizes, an asymptotics FSS 
behavior of the form f3J , 

C max = a c + b c N + ... (9) 

Xmax = a x + b x N + ... (10) 
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Figure 5: Plot of the specific heat C max (circle) and susceptibility Xmax 
(square) versus N. In the Part (a) m — 2 and Part (b) m = 7. 

B hmm = a Bi +bB i N + ... (II) 

We have performed Monte Carlo simulation on directed Barabasi-Albert 
networks with values of connectivity m = 2 and 7. For a given m, we used 
systems of size N = 250, 500, 1000, 2000, 4000, 8000, and 16000 sites. We 
waited 10000 Monte Carlo steps (MCS) to make the system reach the steady 
state, and the time averages were estimated from the next 10000 MCS. In 
our simulations, one MCS is accomplished after all the N spins are updated. 
For all sets of parameters, we have generated 20 distinct networks, and have 
simulated 20 independent runs for each distinct network. 

Results and Discussion 

Our simulations, using the HeatBath algorithm, indicate that the model 
displays a first order phase transition. Fig.l shows the overall behaviour 
of the spontaneous magnetisation. In Fig. 2 we show the dependence of 
the Binder parameter Bi(K) for connectivity m = 2 and 7 on the inverse 
of temperature K and various systems size. Part (a) shows the curves for 
m = 2 of bottom to top of N = 250 to 16000, part (b) the same as part (a) 
for m = 7. The Binder parameter clearly goes to a value which is different 
from 2/3. This is a sufficient condition to characterize a first-order transition. 
In Fig. 3 we plot the Binder parameter B^ versus 1/N for m = 2 (circle) 
and m = 7 (square), and severals systems sizes (N = 1000, 2000, 4000, 8000 
and 16000). We show the scaling of the Binder parameter minima, and again 
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the first order phase transition is verified. The order of the transitions can 
be confirmed by ploting the values of 2/3 — B ijmin again versus 1/N. For 
a second-order transition the curves goes to zero as we increase the system 
size. Here, the quantity 2/3 — -Bi jTO j„ approaches a nonvanishing value in the 
limit of small 1/N as for m = 2 than as m = 7, see Fig. 4. 

As decipted in Figure 5, our results for scaling of the specific heat and 
susceptibility are consistent with equations (9,10). Part a shows m = 2, and 
part b, m = 7. The same occurs with the plot the Fig. 4 for the Binder 
parameter minima, equation (11). In the part a for m = 2 and m = 7 part 
b. 

Conclusion 

In conclusion, we have presented a very simple equilibrium model on 
directed Barabasi- Albert network |2j. Different from the spin 1/2 Ising 
model, in these networks, the spin 1 Ising model presents a the first-order 
phase transition which occurs in model with connectivity m = 2 and m = 7 
here studied. We also verific that occur a phase transition for Potts Model 
for q = 3 and q = 8 on directed Barabasi-Albert network [To] . 
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